Investigation of oscillation frequency and disorder induced dynamic phase transitions 

in a quenched-bond diluted Ising ferromagnet 
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Frequency evolutions of hysteresis loop area and hysteresis tools such as remanence and coercivity 
of a kinetic Ising model in the presence of quenched bond dilution are investigated in detail. The 
kinetic equation describing the time dependence of the magnetization is derived by means of effective- 
field theory with single-site correlations. It is found that the frequency dispersions of hysteresis loop 
area, remanence and coercivity strongly depend on the quenched bond randomness, as well as 
applied field amplitude and oscillation frequency. In addition, the shape of the hysteresis curves 
for a wide variety of Hamiltonian parameters is studied and some interesting behaviors are found. 
Finally, a comparison of our observations with those of recently published studies is represented and 
it is shown that there exists a qualitatively good agreement. 
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I. INTRODUCTION 



When a ferromagnetic material is subject to a time 
dependent oscillating magnetic field, the system may not 
respond to the external magnetic field instantaneously, 
and this situation leads to some interesting phenomena. 
Namely, dynamic phase transitions (DPTs) and dynamic 
hysteresis behavior originate due to a competition be- 
tween time scales of the relaxation time of the system 
and oscillation period of the external magnetic field. On 
the other hand, at high temperatures and for the high 
amplitudes of the periodic magnetic field, the system is 
able to follow the external field with some delay while this 
is not the case for low temperatures and small magnetic 
field amplitudes. This spontaneous symmetry breaking 
indicates the presence of a DPT [l[ which shows itself 
in the dynamic order parameter (DOP) which is defined 
as the time average of the magnetization over a full pe- 
riod of the oscillating field. According to early exper- 
imental studies, dating back to last century, there ex- 
ists an empirical law for the hysteresis loop areal scal- 
ing 0. Up to now, a fairly well numerous theoretical 
studies have been dedicated to DPT and hysteresis prop- 
erties of kinetic Ising models by using various methods 
such as Monte Car lo ( MC) simulations 043) effective- 
field theory (EFT) [11 [13 and mean-field theory (MFT) 
[I. IR [l7l - [T9l . I25l - l29l ] . Apart from these, there exist a de- 
tailed geometrical description [3(| , and also an hysteresis 
criterion based on rate competition between the criti- 
cal curvature and the potential-barrier height [3l| . By 
employing standard Metropolis MC algorithm with peri- 
odic boundary conditions, the pure kinetic Ising model in 
a two dimensional square lattice has been simulated by 
Lo and Pelcovits 0, and they found some evidences of 
a DPT. They also investigated the behavior of dynamic 
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hysteresis for varying Hamiltonian parameters, however 
they did not intend to make any classification for dy- 
namic hysteresis behavior. Afterwards, by making use 
of MC simulations and MFT, Acharyya and Chakrabarti 
[ij have presented a comprehensive investigation of the 
pure kinetic Ising model and they showed that the hys- 
teresis loops are asymmetric in the dynamically ordered 
phase while in a dynamic disordered phase, they are sym- 
metric around the origin in a magnetization versus field 
amplitude plane. In Ref. [5|], the attention has been 
focused on the behavior of dynamic hysteresis and dy- 
namic loop area (DLA) of a pure kinetic Ising model 
by using both MFT and MC calculations. In addition, 
the frequency variations of the coercive field have been 
analyzed below the zero field critical temperatures. It 
has been concluded that the numerical results for the 
dynamic coercive field show a power law frequency vari- 
ation for both cases. However, the coercive field calcu- 
lated by MFT becomes frequency independent in the low 
frequency regime. Differences between the two methods 
(MFT and MC) have been explained by using a kind of 
Landau-type double well free energy mechanism. More- 
over, there exist many previously published studies such 
as the domain wall motion and nucleation theory pll - tzU 
regarding the hysteresis loop mechanism [20] . 

On the experimental picture, it can be said that a 
great many of experimental studies have been devoted to 
DPTs and hysteresis behaviors of Co films on a Cu (001) 
surface [II], [Co(4A°)/Pt(7A°)] multi-layer system with 
strong perpendicular anisotropy [33j , thin polycrystalline 
NigoFe™ films [34[, ultrathin ferromagnetic Fe/Au(001) 
films [HI, epitaxial Fe/GaAs(001) thin films [Hj], epi- 
taxial single ferromagnetic fee NiFe(OOl), fee Co(OOl), 
and fee NiFe/Cu/Co(001) layers J37|, Feo.42Zno.5sF2 [Hf, 
finemet thin films (39|, and epitaxial Fe/GaAs(001) and 
Fc/InAs(001) ultrathin films 0. After some detailed 
experimental investigations, it has been observed that ex- 
perimental non-equilibrium dynamics of considered real 
magnetic systems strongly resemble the dynamic behav- 
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ior predicted from theoretical calculations of a kinetic 
Ising model. It is clear from these works that there ex- 
ists a strong evidence of qualitative consistency between 
theoretical and experimental studies. 

Besides, the lattice models including impurities and 
defects have attracted considerable attention, because 
they are very useful to investigate the behavior of dis- 
ordered systems in nature. For this purpose, more ef- 
fective models have been introduced to examine the in- 
fluences of disorder on the thermal and magnetic prop- 
erties of real magnetic materials. For example, the ef- 
fects of impurities on driving-rate-dependent energy loss 
in a ferromagnet under the time dependent magnetic field 
have been analyzed by Zheng and Li (4lj by using sev- 
eral well defined models within the frameworks of MFT 
and MC, and they found using MFT that, the hysteresis 
loop area is a power law function of the linear driving 
rate as A — Aq oc h 13 , where, Aq, h and /3 are the static 
hysteresis loop area, the linear driving rate and scaling 
exponent of the system, respectively. Very recently, the 
quenched site and bond diluted kinetic Ising models un- 
der the influence of a time dependent oscillating magnetic 
field have been analyzed by making use of EFT |42i. 14311 
on a two dimensional honeycomb lattice. In Rcf. [42j . 
the global phase diagrams including the reentrant phase 
transitions are presented by the authors for site diluted 
kinetic Ising model, and they showed that the coexis- 
tence regions disappear for sufficiently weak dilution of 
lattice sites. Following the same methodology, the au- 
thors have concentrated on the influences of quenched 
bond dilution process on the dynamic behavior of the 
system. After some detailed analysis, it has been found 
that the impurities in a bond diluted kinetic Ising model 
give rise to a number of interesting and unusual phenom- 
ena such as reentrant phenomena and the impurities to 
have a tendency to destruct the first-order transitions 
and the dynamic tricritical point [43| . Furthermore, it 
has also been shown that dynamically ordered phase re- 
gions get expanded with decreasing amplitude which is 
more evident at low frequencies. 

Very recently, frequency dependencies of the dynamic 
phase boundary (DPB) and magnetic hysteresis of Ising 
model in an oscillating magnetic field have been exam- 
ined by Punya et al. [19| . Based on MFT, they have 
represented the phase diagram in a magnetic field ver- 
sus temperature plane for varying frequencies, in order 
to show the influence of the oscillation frequency of the 
external field on the hysteretic response of the system. 
According to their calculations, it has been found that 
the frequency dispersion of the DLA, the remanence and 
coercivity can be categorized into three distinct types 
for a fixed temperature. Although the detailed investi- 
gations and classifications of the hysteretic behaviors of 
pure kinetic Ising model have been done within the MFT 
|19(, as far as we know, the dynamic hysteretic behavior 
within the EFT has not been studied for disordered ki- 
netic Ising models. Actually, the main motivation for our 
work comes from the recent analysis of Ref. [l9j for the 



investigation of hysteresis behavior and dynamic phase 
diagrams of kinetic Ising model. Therefore in this work, 
we intend to probe the effects of the quenched bond dilu- 
tion process on the dynamic hysteresis behavior of kinetic 
Ising model in the presence of a time-dependent oscillat- 
ing external magnetic field by using the EFT with corre- 
lations based on the exact Van der Waerden identity for 
a spin-I/2 system. Present method is quite superior to 
MFT, since thermal fluctuations are partly considered in 
EFT. The outline of the paper can be summarized as fol- 
lows: The dynamic equation of motion, DOP and DLA 
of the quenched bond diluted kinetic Ising model have 
been introduced in the next section. The numerical re- 
sults and related discussions are given in Section Unl and 
finally Section llVI contains our conclusions. 



II. FORMULATION 

The kinetic Ising model is given by the time dependent 
Hamiltonian 



* = - E J " 



QZ QZ 



(1) 



<ij> 



where the spin variables = ±1 are defined on a hon- 
eycomb lattice (q = 3) and the first sum in Eq. JT]) is 
over the nearest neighbor pairs of spins. We assume that 
the nearest-neighbor interactions are randomly diluted 
on the lattice according to the probability distribution 
function 



P(J iS ) = p8{Jij - J) + (1 - p)S{Jij), 



(2) 



where p denotes the concentration of active bonds. The 
term H(t) in Eq. ([1]) is a time dependent external mag- 
netic field which is defined as 



H(t) = hcos(ujt), 



(3) 



where h and w = 27r/ represent the amplitude and the 
angular frequency of the oscillating field, respectively. If 
the system evolves according to a Glauber-type stochas- 
tic process [44| at a rate of 1/r which represents the tran- 
sitions per unit time then the dynamic equation of motion 
can be obtained as follows 



r^ = -(Sf) + (tanh 



Ei + H{t) 



(4) 



where Ei = Y^j JijSj is the local field acting on the lat- 
tice site i, and ks and T denote the Boltzmann constant 
and temperature, respectively. If we apply the differen- 
tial operator technique (4|| Hg] in Eq. Q by taking into 
account the random configurational averages we get 



dm 
~dt 




/(z)|x=o, (5) 



where A^j — cosh(Jj, V), Bij — sinh(Jj,-V), and m = 
((Sf)) r represents the average magnetization. V = d/dx 
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is a differential operator, and the inner (...) and the outer 
(...),. brackets represent the thermal and configurational 
averages, respectively. When the right-hand side of Eq. 
([5]) is expanded, the multispin correlation functions ap- 
pear. The simplest approximation, and one of the most 
frequently adopted is to decouple these correlations ac- 
cording to 

((S? S?...Sf » p = «S?» r «5?» r ... «S?» r , (6) 

for i j =/=... ^ I [47] . If we expand the right-hand side 
of Eq. (0) with the help of Eq. © then we obtain the 
following dynamic effective-field equation of motion for 
the magnetization of the quenched bond diluted kinetic 
Ising model 

dm 9=3 

i=0 

where the coefficients Xi can be easily calculated by 
employing the mathematical relation exp(aV)/(x) = 
f(x + a). Eq. ([7]) describes the non-equilibrium behavior 
of the system in the EFT formalism. Moreover, the time 
dependence of magnetization can be one of two types ac- 
cording to whether they obey the following property or 
not 

m(f) = -m(t + 7r/w). (8) 

A solution satisfying Eq. ((SJ) is called symmetric solu- 
tion and it corresponds to a paramagnetic or disordered 
phase. In this type of solution, the time dependent mag- 
netization oscillates around zero value. The second type 
of solutions which does not satisfy Eq. © is called non- 
symmetric solution which corresponds to ferromagnetic 
(i.e. ordered) phase where the time dependent magneti- 
zation oscillates around a non-zero value. 

The time averaged magnetization over a full cycle of 
the oscillating magnetic field acts as DOP which is de- 
fined as follows [17| 

Q = ^ / m(t)dt, (9) 

where m(t) is a stable and periodic function. On the 
other hand, the DLA corresponding to energy loss due 
to the hysteresis is defined as [f| 

A = - I m[t)dH = hoj j> m(t) sin(wi)dt. (10) 

The behavior of DOP and DLA as a function of the tem- 
perature for selected Hamiltonian parameters can be ob- 
tained by solving Eqs. jTJ), © and (ITU1) numerically. 
The remanent magnetization and coercivity can be cal- 
culated by benefiting from hysteresis loop (m(t) — H{t) 
curve) in order to understand and clarify the behavior of 
the quenched bond diluted system. 



III. RESULTS AND DISCUSSION 

In this section, in order to understand the dynamic 
evolution of the magnetic system in detail, we will fo- 
cus our attention on the hysteretic response and DPT 
properties of the quenched-bond diluted kinetic Ising 
model under a time-dependent applied magnetic field. 
It is well known that depending on the value of the 
applied field frequency, the magnetic system can exist 
in a dynamically ordered or disordered phase. From 
this point of view, we give the DPB in a (kgT c / J — lu) 
plane with four selected amplitudes of the applied field 
(h = 0.25,0.5,0.75, 1.0), and for some selected values of 
the active bond concentration in Figs. [IJa-d), in order 
to clarify the influence of applied field frequency on the 
DPT properties. At first sight, one can easily observe 
from Figs. Ufa-d) that the dynamically ordered phases 
get expanded with increasing applied field frequency for 
all considered active bond concentrations. Actually, this 
is an expected result since as the frequency increases then 
the time-dependent magnetization can not follow applied 
magnetic field immediately, hence the dynamically or- 
dered phase region moves upward with the increasing 
frequency. We can also mention that aforementioned sit- 
uations are also reported in previously published works 

Qui nam si a 

On the other hand, as the value of p decreases, start- 
ing from pure case (i.e. p = 1.0), the magnetic system 
exhibits unusual and interesting behaviors. For instance, 
dynamically disordered phases occur in the low temper- 
ature regions and also reentrant behavior emerges. An- 
other important consequence of decreasing p is that the 
dynamically ordered phases get narrower for all values 
of the considered applied field amplitudes. These ob- 
servations can also be explained by the following mech- 
anism: The energy contribution which comes from the 
spin-spin interaction gets smaller when the value of p 
is decreased. So, the system can undergo a DPT at a 
lower critical temperature, since the energy contribution 
originating from the temperature and (or) magnetic field 
overcomes the ferromagnetic spin-spin interactions. As a 
result of this mechanism, the dynamically ordered phase 
regions in (ksT c / J — u) plane get narrower. As discussed 
above, the existence or absence of the reentrant behavior 
strongly depends on the applied field amplitude, as well 
as oscillation frequency. As an example, let us consider 
the curves with p = 0.85 to show the influences of the 
applied field amplitude on the DPB, for selected values 
of h. As seen in Figs. Ha-d), f° r h = 0.25 and 0.5, the 
magnetic system exhibits reentrant behavior whereas this 
behavior disappears with further increase in h, such as 
for h = 0.75 and 1.0. 

It is known that applied field frequencies naturally af- 
fect the dynamic evolution of a magnetic system. In order 
to show the influences of the applied field frequencies on 
the system, the temperature dependencies of the DOP 
(Q) and DLA (A) are depicted in Figs. HJa-d) for values 
of p = 0.85 and h = 0.25 corresponding to Fig. [TJa). 
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FIG. 1: Dynamic phase boundaries of the quenched bond diluted kinetic Ising system in a (fcsTc/J — w) plane with some 
selected values of amplitude of the applied magnetic field, (a) h = 0.25, (b) h = 0.50, (c) h = 0.75 and (d) h = 1.0. The 
numbers accompanying each curve denote the value of active bond concentrations. 



In accordance with this purpose, we separated the fre- 
quency space into four regions: The order of the mag- 
netic system is dynamically paramagnetic in first region 
for u> = 0.01. One can see from Fig. [2ja) that the DLA 
gradually increases and exhibits a local maxima with a 
Schottky-like rounded hump and then tends to diminish 
with increasing temperature. On the other side, it is clear 
from Fig. [2f b) that there exists two successive second or- 
der DPTs for a selected frequency value of u> = 0.1, and 
this type of treatment corresponds to second region. The 
first DPT is from paramagnetic to ferromagnetic phase 
while the second DPT is from ferromagnetic to paramag- 
netic phase. For the same Hamiltonian parameters, the 
DLA exhibits two local maxima and two whiskers and the 
positions of the whiskers may correspond to DPT points. 
In Fig. [He), we represent the thermal variation of the 
DOP and DLA curves for u> = 0.5 which indicates the 
third region. It can be easily seen from Fig. [2tc) that 
only one DPT emerges in magnetic system for u> = 0.5. 
DOP curve for ui = 0.5 increases form its ground state 
value with increasing temperature which originates from 
thermal agitations and it decreases with further increas- 
ing temperatures and then falls to zero at a second order 
DPT temperature. Similarly, with the increasing temper- 
ature, starting from zero, the DLA decreases and shows 
a local minima at the dynamically ordered phase region 
and then it exhibits a sharp peak near the DPT point. In 
the following analysis, let us consider the higher value of 
the applied field frequency such as u = 10.0 which refers 
to fourth region. As we can see from Fig. [U[d) that the 



magnetic system exhibits a second order DPT. However, 
we should also mention that the value of phase transition 
temperature observed for to = 10.0 is slightly higher than 
that of to = 0.5. Actually, this is an expected result be- 
cause an increase in the applied field frequency makes the 
occurrence of phase transition difficult, hence the dynam- 
ically ordered phase gets expanded with the increasing 
applied field frequency. For given Hamiltonian parame- 
ters, the DLA shows a peak above the second order DPT 
point. We should note that the existence of a peak on the 
DLA near DPT point has also been indicated in previous 
investigations for various models [f| [H, [H, EE SH • 

We give the frequency dispersion curves of the DLA 
in Figs. [3fa-d) for four selected p values with a variety 
of h values. We should note that the reduced tempera- 
ture in this study is fixed as ksT/J = 1.0 which is lower 
than the transition temperature T c of a pure system with 
EFT for q = 3 in the absence of external field. This selec- 
tion of the temperature allows the system to undergo an 
amplitude, frequency or disorder induced dynamic phase 
transition. Fig[3][a) is plotted for p — 1.0 where all fer- 
romagnetic spin-spin exchange interactions (J > 0) are 
active. At first sight, one can clearly observe in Fig. [SJa) 
that DLA curves can exhibit two types of characteris- 
tic behavior. The first one corresponds to a dynamically 
ordered phase where the DLA increases as the applied 
field frequency increases and then exhibits a character- 
istic rounded peak (or local maxima) for the relatively 
low field amplitude values such as h = 0.1,0.2,0.3,0.4 
and 0.5. These curves are identical to type-I curves of 
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Ref. [19]. As is expected, the DLA is relatively close 
to zero at the low and high applied field frequency re- 
gions in the dynamically ordered phases. On the other 
hand, for h = 0.6, 0.7, 0.8, 0.9 and 1.0 values, DLA curves 
start from a nonzero value and increases with increasing 
frequency values then exhibit a finite cusp at a second 
order phase transition temperature at which the shape 
of hysteresis loop changes from symmetric to asymmet- 
ric and also a rounded hump in high frequency regime 
which corresponds to second type of characteristic be- 
havior. These curves correspond to type-II curves in Ref. 
[l9j . Finite cusps are found to be discontinuous for low 
h and weak disorder. The rounded humps observed in 
DLA curves in high frequency regime disappears for in- 
creasing h and (or) decreasing p values. For instance, for 
h = 0.6 with p — 1.0, the DLA exhibits a discontinu- 
ous cusp and a rounded distinct hump while it displays 
only a continuous cusp for ft, = 1.0 and p = 1.0. We can 
mention that the characteristic cusp of (h,p) = (0.6, 1.0) 
curve arises from the frequency evolution of symmetric 
hysteresis loop, whereas the rounded hump of this curve 
originates from the frequency evolution of asymmetric 
hysteresis loop. On the other side, it is known that the 
phase transition temperature decreases with increasing 
amount of disorder. The effects of the aforementioned 
situation on the frequency dispersion can be seen in Figs. 
[3Kb-d) which have been plotted for p = 0.9, 0.8, and 0.7. 
One can observe in Figs. [3][b-d) that as the disorder ef- 
fects become dominant in the system then the representa- 
tive type-I curves tend to evolve into type-II, and type-II 
curves have a tendency to turn into type-Ill curves de- 
fined in Ref. [l9| which represent the paramagnetic DLA 



dispersion curves. In contrast to type-I curves, type-Ill 
curves have a finite residual value in the limit of uj — > 0. 
In addition, the type-Ill curves are characterized by a 
rounded hump, as type-I curves. Based on these obser- 
vations, we see that DLA dispersion curves clearly reveal 
the signs of the oscillation frequency and disorder induced 
dynamic phase transitions in the present system. 

In Figs. EJa-d), we present the frequency dependence 
of remanent magnetization curves corresponding to DLA 
curves in Fig. [3] Here we select four typical values of dis- 
order parameter p with a variety of field amplitude h val- 
ues. Remanent magnetization is defined as the residual 
magnetization of a ferromagnetic material after an ex- 
ternal magnetic field is removed. Hence the remanence 
data can be easily extracted from calculated hysteresis 
loops. At first glance, it can be easily seen that depend- 
ing on amplitude and frequency of the external magnetic 
field, as well as the amount of disorder, remanence curves 
can exhibit three distinct characteristics corresponding to 
type-I, type-II and type-Ill profiles of DLA curves. Cor- 
responding to typc-I DLA curves, either branches of re- 
manent magnetization curves are positively valued which 
means that the system is ordered ferromagnetically in the 
whole range of oscillation frequency u>. Increasing the 
amount of disorder or amplitude of the external field, rep- 
resentative type-I remanence curves may evolve to type- 
II and type-Ill curves in order. Type-II remanence curves 
exhibit equally valued branches with opposite signs in 
the low frequency region indicating that the system is in 
a dynamically paramagnetic state. As the frequency in- 
creases then the negative branch of a representative type- 
II remanence curve changes its sign at a critical value of 
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FIG. 3: The frequency dispersion of the hysteresis loop area curves with considered applied field amplitudes for selected values 
of active bond concentrations (a) p = 1.0, (b) p = 0.9, (c) p — 0.8 and (d) p — 0.7. 



frequency, hence we understand that the corresponding 
hysteresis loop loses its symmetry and becomes dynam- 
ically ordered. Hence this behavior can be considered 
as a frequency induced phase transition. For sufficiently 
high amplitudes and (or) for a sufficiently large amount 
of disorder, the observed type-II curves evolve to type- 
Ill characteristics. In this case, either branches of re- 
manence curves are symmetric with opposite signs which 
indicates that the system is dynamically disordered for 
the whole range of oscillation frequency spectrum. Re- 
manent magnetization is a vital property for materials 
used in magnetic information storage devices. Namely, 
in case of some fixed values of external field amplitude 
and oscillation frequency, as well as an amount of exist- 
ing disorder, for a device made of a material with higher 
remanence property there will be less chance that the 
information stored in the media will be erased due to 
some external sources [4!|. From this point of view, it 
is important to note that according to Fig. 0|d), there 
exists an optimum field frequency for a disordered ma- 
terial under an oscillating magnetic field at which the 
type-Ill remanence (i.e. the distance between positive 
and negative remanent magnetizations) exhibit a maxi- 
mum value, and as the field amplitude increases then the 
optimum frequency value shifts towards higher frequency 
values. 

As another characteristic feature of hysteresis loops, 
we have examined the coercivity of the system which is 
defined as the minimum required reversal magnetic field, 
in order to reduce the remanent magnetization to zero. 
For storage media, a material with a high coercivity is 



desired to preserve the stored data from unwanted exter- 
nal factors. We represent the variation of the coercive 
field curves as a function of oscillation frequency in Figs. 
[5Ja-d) for some selected values of disorder parameter p 
and external field amplitude h. For the DLA curves with 
type-I characteristics, there is no coercivity for any field 
amplitude, oscillation frequency, or disorder, because all 
hysteresis curves corresponding to type-I behavior have 
asymmetric shape. This means that the instantaneous 
magnetization can not follow the external magnetic field. 
For type-II characteristics however, we see that at low 
frequencies, absolute values of either branches in a coer- 
civity curve of the system increase for a while then tend to 
decrease with further increasing frequencies. After a crit- 
ical frequency value, these branches lose their symmetric 
shape then coercive fields disappear and a dynamic phase 
transition from a dynamic paramagnetic phase to a fer- 
romagnetic one takes place. We note that nonsymmetric 
coercivity data have not been shown in Fig. [5] We again 
see that an oscillation frequency induced dynamic phase 
transition takes place in representative type-II curves. If 
the magnetic field amplitude is large enough or if there 
exists a sufficiently strong disorder (see Figs. [SJ; and 
d), coercivity curves may evolve from type-II to type-Ill 
characteristics. Since the type-Ill hysteresis curves are 
symmetric, coercive field can be observed for the whole 
range of oscillation frequency values. However, for type- 
Ill behavior, coercivity is rather small at low frequencies. 
As the field frequency increases then coercivity of the sys- 
tem increases and in the limit of to — > oo it saturates to 
its maximum value. The maximum value of coercivity 
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FIG. 4: Frequency dispersion of remanent magnetization for a variety of applied field amplitudes with selected active bond 
concentrations (a) p = 1.0, (b) p = 0.9, (c) p = 0.8 and (d) p = 0.7. 



is noting but just the applied amplitude of the external 
magnetic field for sufficiently high frequencies. Hence, we 
can conclude that the greater amplitude of the external 
magnetic field, we get more enhanced coercivity. 

The two fundamental tools, namely coercivity and re- 
manence properties allow us to make predictions on the 
shape of the hysteresis loops. Effect of the field fre- 
quency on the shape of typical hysteresis loops corre- 
sponding to Figs. |3l |H and [5] are depicted in Fig. [6] for 
a wide variety of oscillation frequencies within the range 
0.01 < ui < 10.0 with given values of p and h. In Fig. 
Ella), we present the hysteresis loop curves correspond- 
ing to type-I characteristics of DLA curves. As seen in 
Fig. [^a), hysteresis curves exhibit positive valued nar- 
row loops at low frequencies then evolve into a horizon- 
tal line with increasing frequencies. Type-II curves are 
shown in Fig. [(Jb), where the low frequency paramag- 
netic loops turn into high frequency ferromagnetic loops 
due to a frequency induced dynamic symmetry loss in 
the system. Moreover, according to Fig. [He), symmetric 
type-Ill hysteresis curves keep their dynamic paramag- 
netic state and exhibit remanence and coercivity with 
branches having both positive and negative values. We 
can conclude from Fig. |6jc) that low frequency thin 
and iS-shaped hysteresis loops evolves into broad square- 
like shapes then turn into elliptical loops for sufficiently 
high frequencies. If we increase the frequency further, 
these elliptical loops become horizontal lines. Related 
to the frequency evolutions of the hysteresis curves, we 
should also note that a qualitative consistency is found 
between our numerical findings and experimental obser- 



vations m musl- 
in the following analysis, let us focus on the influence 
of disorder on the frequency dispersion of hysteretic prop- 
erties of the system. In Figs. Ufa) and[7Jb), we repre- 
sent variation of coercivity and remanence curves as a 
function of oscillation frequency for some selected val- 
ues of disorder parameter p with a given field amplitude 
h = 0.5. From Fig. [7[a), we see that as the concentra- 
tion of diluted bonds increases then the zero frequency 
branches of coercivity curves vanish as u> — » 0. In type-II 
curves, critical frequency at which the system exhibit a 
dynamic phase transition increases with increasing dis- 
order. This frequency value can be regarded as a cut-off 
frequency at which the symmetric loops turn into asym- 
metric loops. Hence, the non-symmetric coercivity data 
has not been considered in Fig. [TJa). Moreover, coerciv- 
ity curves exhibiting type-II characteristics may evolve 
into type-Ill behavior, if the concentration of diluted 
bonds exceeds a critical value (i.e. percolation thresh- 
old). Hence, we understand that the system exhibits a 
disorder induced dynamic phase transition. In Fig. EJb), 
we illustrate the dependence of frequency dispersion of re- 
manence curves on the concentration of disordered bonds 
with a fixed external field amplitude h = 0.5 correspond- 
ing to Fig. Efa). Type-I curves in Fig. (7{b) exhibit a pos- 
itive remanence which indicates that the system stays in 
a dynamic ferromagnetic order for the whole range of fre- 
quency spectrum. By increasing disorder, type-I curves 
turn into type-II characteristics where we can observe 
that either branches of remanence curves tend to lose 
their symmetric shape with increasing frequency which 
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FIG. 5: Frequency dispersion of the coercivity for a variety of applied field amplitudes with selected values of active bond 
concentrations (a) p = 1.0, (b) p = 0.9, (c) p = 0.8 and (d) p = 0.7. 



signals the appearance of a frequency induced dynamic 
phase transition. In a similar manner, we observe in Fig. 
[7](b) that for sufficiently low concentration of magnetic 
bonds, the remanent magnetization curves may exhibit 
type-Ill characteristics where the system always exhibit 
a dynamically disordered phase for the whole range of 
oscillation frequency spectrum. 



As a final investigation, we discuss the evolution of 
some typical hysteresis loops corresponding to remanence 
and coercivity properties depicted in Fig. [7J In Fig. [5] 
we fix the amplitude of the external magnetic field as 
h = 0.5. Fig. [Sfa) corresponds to low frequency situa- 
tion (cj = 0.01). From Fig. [5£i, we see that a disorder 
induced dynamic phase transition takes place in the sys- 
tem with increasing disorder. Hysteresis loops exhibit a 
positive valued horizontal line in the pure case whereas as 
the amount of disorder increases then the loop area gets 
narrower and we observe S-shaped thin loops. On the 
other hand, the situation is depicted for relatively high 
frequency values such as w = 1.0 in Fig. [8Kb). In this 
case, asymmetric loops corresponding to ferromagnetic 
order become elliptical symmetric loops in the dynamic 
paramagnetic phase with increasing disorder. Hence, we 
can conclude that influence of disorder on the shape of 
hysteresis loops strictly depends on the magnitude of the 
applied field frequency. Aforementioned results are also 
consistent with the observations illustrated in Fig. [7] 



IV. CONCLUSIONS 

In conclusion, by using EFT with correlations, we 
have investigated frequency dispersion features of dy- 
namic loop area, coercivity and remanent magnetization 
of a kinetic Ising model in the presence of quenched- 
diluted bonds. First of all, in order to gain some gen- 
eral insight regarding the influence of bond disorder, ex- 
ternal field amplitude, and oscillation frequency on the 
dynamic phase transition properties of the system we 
have demonstrated the phase diagrams of our model in 
a (ksTcj J — lo) plane for weak (h = 0.25), intermediate 
(h = 0.5,0.75) and large (h = 1.0) values of the external 
field amplitude with some selected values of disorder pa- 
rameter p. We have found that critical frontiers between 
symmetric and asymmetric phases in (fc,BT c / J — u) plane 
are strongly dependent on amplitude and frequency of os- 
cillating field, as well as amount of disorder. 

EFT is superior to conventional MFT predictions ow- 
ing to the fact that the single spin correlations, in 
other words thermal fluctuations are partially considered. 
Hence, we take the recently published work performed by 
Ref . [19( one step forward by considering single-site corre- 
lations between different lattice sites, and by considering 
the presence of disordered bonds. However, it is worth to 
note that the approximation method used in the present 
work, does not take into account the local structure of the 
lattice. In other words, EFT cannot distinguish between 
a two dimensional triangular lattice and a three dimen- 
sional simple cubic lattice, since both have a coordination 
number q = 6. In this context, the results obtained in 
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FIG. 6: Hysteresis loops of the quenched bond diluted kinetic 
Ising system for several values of applied field frequencies and 
amplitudes with some selected values of active bond concen- 
trations. 



the present paper depend on the coordination number of 
the system. 

The results obtained in the present paper qualitatively 
support the recent findings based on MFT [l9[ for the 
pure case. Namely, we have observed that frequency dis- 
persion of dynamic loop area (DLA) can be classified 
into three distinct types. In type-I behavior, the system 
always exhibits asymmetric hysteresis loops at any os- 
cillation frequency for small field amplitudes and weak 
amount of disorder whereas type-II curves are charac- 
terized by a frequency induced phase transition where 
the low frequency symmetric loops become asymmetric 
loops with increasing oscillation frequency. On the other 
hand, type-Ill behavior is essentially represented by sym- 
metric hysteresis loops for the whole range of frequency 
spectrum when the external field amplitude is sufficiently 
large and (or) if there exists sufficiently strong dilution 
of ferromagnetic nearest-neighbor bonds. Therefore we 
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FIG. 7: Variation of the coercive field (a) and the remanent 
magnetization (b) as a function of the applied field frequency 
for a field amplitude value of h = 0.5 with some selected 
values of disorder parameter p. 
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FIG. 8: Hysteresis curves of the quenched bond diluted ki- 
netic Ising system for several values of active bond concen- 
trations with selected values of applied field frequencies and 
amplitude. 
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can conclude that type-II curves exhibit frequency in- 
duced dynamic phase transitions whereas evolution of 
DLA curves from type-I to type-Ill characteristics is due 
to a disorder induced effect or as a consequence of large 
amplitudes of the external magnetic field. 

Aforementioned observations have also been supported 
by detailed investigation of frequency dispersions of co- 
ercivity and remanent magnetizations. From the techno- 
logical point of view, in magnetic information storage de- 
vices a material with high coercivity and large remanence 
is desired to preserve the stored data from unwanted ex- 
ternal factors. In this regard we have found that there 
exists an optimum field frequency for a disordered ma- 
terial under an oscillating magnetic field at which the 
type-Ill remanence (i.e. the distance between positive 
and negative remanent magnetizations corresponding to 
type-Ill) exhibit a maximum value, and as the field am- 
plitude increases then the optimum frequency value shifts 
towards higher frequency values. Finally, we have ob- 



served that influence of disorder on the shape of hystere- 
sis loops strictly depends on the magnitude of the applied 
field frequency. 

As a conclusion, we hope that the results obtained in 
this work would shed light on the further investigations 
of the dynamic nature of the critical phenomena in dis- 
ordered systems and would be beneficial from both the- 
oretical and experimental points of view. 
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